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Abstract 



We use the experimental limit on the interference of Ml and E2 multipoles 
in the 7-decay of 57 Fe to bound the time-reversal-violating parity-conserving 
pNN vertex. Our approach is a large-basis shell-model calculation of the 
interference. We find an upper limit on the parameter g p , the relative strength 
of the T- violating pNN vertex, of close to 10 -2 , a value similar to the best 
limits from other kinds of experiments. 

PACS: 21.60.Cs, 23.20.Gq, 24.80.+y 
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For many years it has been difficult to compare the quality of limits on time-reversal- 
violating parity-conserving (TVPC) interactions coming from different low-energy experi- 
ments. The experiments typically limit observables unique to themselves, and before com- 
parisons can be made, these limits must be translated into a common TVPC quantity. It 
turns out that a convenient measure of nuclear TVPC interactions is the dimensionless ratio, 
often called g p [fj]], of the TVPC p- meson - nucleon coupling to the normal strong coupling 
g p . Among the other mesons only those with axial- vector couplings can transmit TVPC in- 
teractions between nucleons via a single exchange || , and they are significantly heavier than 
the p and consequently less effective in nuclei. It is therefore reasonable to treat all TVPC 
nucleon-nucleon interactions as arising from p exchange, and to use g p to parameterize their 
strength. 

Experimental upper limits on several quantities, including the electric dipole moments 
of the neutron and of 199 Hg |]], and a correlation in the scattering of polarized neutrons 
from aligned 165 Ho ||, have been translated into limits on g p , constraining it to be less 
than about 10 -2 . A number of other experiments, looking e.g. for the violation of detailed 
balance J3J, remain to be similarly interpreted. In this paper we report an examination of 
a 1977 experiment that searched for interference between Ml and E2 radiation in the 
7-decay of the first 5/2~ state in 57 Fe to the first 3/2 _ state. (Neither is the ground state; 
the two have excitation energies of 137 keV and 14 keV.) Our approach was to diagonalize 
the strong nuclear hamiltonian in the shell-model, and then treat the TVPC p-exchange 
interaction as a perturbation that causes the interference by mixing higher-lying states into 
the two involved in the transition. Ref. || employed this method to constrain the TVPC 
coupling of the A\ meson to the nucleon from the same experiment, but used what we argue 
is too small a model space. In addition, the lighter and more commonly considered p meson 
was neglected completely. 

The M1-E2 interference that signals T violation can be expressed in terms of sin^PJ the 
imaginary, part of the multipole mixing ratio S M\, which is defined as M 
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In Ref. || the upper limit on | sin rj | was expressed in the form of a measured value that 
included zero within experimental accuracy: 

I sin 77 1 = (3.1 ±6.5) x 1(T 4 . (2) 

The contributions to 77 can be written as 

V = e E2 - e M i + £ £S2,£mi,£<1 (3) 

where the last term £ is represents effects of final state interactions, which have been shown 
to be smaller than the upper limit in Eq. (0). 



1 sinr/ is directly proportional to the experimental correlation (J-qxE)(J-q)(J-E), where J is the 
quantization axis of the initial nucleus, E is the photon electric field vector, and q is the photon 
direction ||. 
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In first-order perturbation theory, the difference between the two e's is 
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With nucleons represented by i, j, the two-body p-exchange potential has the form 

v p = EK ^ XT ^3 (5) 



-m ri 



where r^- = — 1 = r^- x -(p, — p^), = 3.70 n.m. is the isovector nucleon magnetic 
moment, M is the nucleon mass, g p = 2.79 is the normal strong pNN coupling, and g p is 
the dimensionless ratio of the TVPC coupling to g p that we are trying to constrain. After 
choosing a model space and interaction (and a reasonable prescription for treating short- 
range correlations fHH), we can use this formalism in a shell-model calculation to translate 
the experimental limit on sin t] to a limit on g p . 

The issues surrounding the calculation are more complicated than they initially appear, 
however. To evaluate the phases in Eq. one needs, in principle, the wave functions and 
energies of all 3/2~ and 5/2~ states in 57 Fe. To obtain them, one ought to diagonalize the 
best available nuclear hamiltonian for 17 valence nucleons moving freely in the pf shell. 
Such a space has an m-scheme dimension of ~ 4.5 x 10 8 . At the other extreme is a minimal 
model space, based on the well-established shell closure at iV or Z = 28, consisting of 3 
valence neutrons in the (2p 3 / 2 , 1/5/2; ^Pi/2) shells and the remaining 14 nucleons in the I/7/2 
subshell. This "small space" is the one used in Ref. and contains few enough states 
to allow direct diagonalization of any Hamiltonian. Unfortunately this space artificially 
restricts the Ml strength from any given state because it doesn't allow the important I/7/2- 
I/5/2 spin-flip transition. Consequently, in the calculations described here we used a "large 
space", constructed by allowing a single proton or neutron to move out of the I/7/2 shell 
into any one of the other subshells. The large space contains 23604 m-scheme states, forcing 
an approximate diagonalization. 

To obtain the approximate wave functions we used the Lanczos algorithm as implemented 
in the shell-model code CRUNCHER |nj and its auxiliary codes, with slight modifications 
to accommodate the imaginary two-body matrix elements of the interaction in Eq. (^j). 
Since it was not practical to calculate all J n = | and | wave functions (there are 2052 
and 2755 of these, respectively), we adopted a procedure expounded in Ref. JL2] to obtain 
Gamow- Teller strength functions. We first used the Lanczos algorithm to obtain the lowest 
I and I states in 57 Fe to high precision. Next we created a "collective" E2 or Ml state 
by acting on the parent state with the relevant operator. We then used the collective 
state as the initial basis vector for an approximate Lanczos-based diagonalization of higher- 
lying states, yielding pseudo-eigenvectors (PSEVs), which approximate the true states. We 
typically performed about 100 Lanczos iterations, resulting in about 100 PSEVs for each of 
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the Jj's. In the Lanczos approach the lowest (and highest) several PSEVs are quite accurate 
representations of the corresponding eigenstates, while at intermediate energies the PSEV's 
converge more slowly, and after pa 100 iterations each still has contributions from tens to 
hundreds of actual eigenstates. It is easy to see, however, that all of the strength is contained 
in these PSEVs, which we used for the states nJi and nJf in Eq. (Q). The wave functions 
used in each of the four terms in Eq. were slightly different since they originated from 
different initial collective states. 

So far we have not mentioned our choice of interaction. There are several effective 
interactions on the market, but (unfortunately) we did not know which was the best in 
this space. We were able to test the sensitivity of our results to the choice of Hamiltonian, 
however, and so used 3 different p /-shell interactions: the FPVH interaction of [ 13fl , the 
TBLC8 interaction of ||I4" |, and the FPBPN interaction (the FPD6 interaction of |15| with 
the single particle energies modified to fit 56 Ni |16|)- Each of these interactions reproduced 
the energy spectrum of low lying states in 57 Fe reasonably well. The spread in the calculated 
values of the phase £e2~ £ mi with these interactions provided a rough measure of theoretical 
uncertainty. 

The last component of the calculation was the choice of effective E2 and Ml operators 
for each force. The matrix element (||| E2 or Ml |||) normalizes each term in Eq. (^). Since 
the Ml matrix element (in the denominator in the second and fourth term in Eq. (^)) is very 
small, it is particularly important, and we chose effective (/-values for the Ml operator in 
order to reproduce it accurately. Our prescription was to fix all of the Ml ^-values, except 
for the isoscalar spin piece {g s is), at their free nucleon values. For each interaction we then 
chose g s IS to give the correct matrix element for the first transition. The sign of the matrix 
element is not known, so we chose it consistently amongst the forces to obtain the most 
reasonable values for the set of g s IS 's. 

For the E2 operator a similar procedure gave unrealistic values for the effective charges 
e p and e n ; we therefore adopted the "canonical" values e p = 1.5e and e n = 0.5e for all 
of the interactions. These values result in reasonable agreement with the first E2 matrix 
element, especially for the FPBPN force. In addition, the final phase e^ 2 is only weakly 
dependent on the choice of the E2 effective charges. Table I summarizes the E2 and Ml 
matrix elements and total strengths for the few lowest states in both the large and small 
spaces. (The TBLC8 force shares a common heritage with the FPVH force and, since the 
results are similar, we omit TBLC8 from the tables). In the large space the total strength 
for both multipoles is relatively insensitive to the force chosen. However, the Ml strength is 
about a factor of 10 larger than in the small space, dramatically illustrating the importance 
of including the I/5/2 level. 

How much did the non-convergence of the intermediate PSEVs affect the results? The 
answer is very little for the E2 part of the phase, because the strength is concentrated at 
low energies and the energy denominator in Eq. (^) enhances the contribution of the low 
lying converged states and reduces the effects of the higher lying states. In Figure 1 we show 
the distribution of E2 strength for the | — > n| transitions (dashed line). It is completely 
dominated by transitions among the converged states. A similar result holds for the E2 in 
the I — > n| direction. Though the effects of the energy denominator are also at work in 
Ml piece of the phase, the distribution of Ml strength complicates matters. In Figure 2 we 
show the total Ml strength for the § — > n| transitions (dotted line); a broad resonance 
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is visible at ~ 12 MeV. Although this is the region where the PSEVs are unconverged, the 
Ml part of the phase nonetheless seems to be represented reasonably well. We make this 
statement after varying the number of Lanczos iterations and hence the number of PSEVs 
(converged and unconverged) to see if the phase changed appreciably as the approximations 
became more accurate. The size of the dependence is illustrated in Table II, where the E2 
and Ml parts of the phase t] (with g p — 1) are listed for several numbers of iterations and 
for two different interactions. The E2 phases show essentially no dependence on the number 
of iterations (as implied above) and the Ml phases are not affected dramatically, indicating 
that the true result is not far from our best approximation. 

The results in Table II allow us to constrain the parameter g p and estimate the uncer- 
tainty. The FPVH and the FPBPN forces give very similar results for each piece of the 
phase and the final phases are very close. The TBCL8 interaction gives a similar result, 
£E2 — £ mi — —24.2 x 1CT 3 . The lack of dependence on the interaction suggests that the 
uncertainty in the results is not large. Table II also suggests that the phase is insensitive to 
the size of the model space, but this turns out to be a coincidence. In the small space, all of 
the Ml piece of the phase lies at very low excitation energy, mirroring the initial downward 
peak at 2-3 MeV in Figure 2. But the fall in the phase from 3-10 MeV and the subsequent 
rise due to the Ml resonance are not present in the small space and so the agreement on 
the final value of Emi between the two model spaces is accidental. 

The entries in Table I were evaluated with g p — 1. Neglecting theoretical error, which we 
have argued should be fairly small, and averaging the results from the FPVH and FPBPN 
forces in the large space, we conclude that \ee2 — £Mi\/g P = 16.4 x 10~ 3 . The experimental 
value for | sin^|, Eq. (0), then implies that 

|^| = (2±4)xl0" 2 . (6) 

This number is comparable to the best limits from other experiments. Limits on electric 
dipole moments, for example, correspond to \g p \ ^ 10~ 2 , and the new data on neutron- 
holmium [17| scattering yields \g p \ = (2.3 ± 2.1)10~ 2 . Perhaps coincident ally, all these very 



different experiments give roughly the same limit. It has been suggested ||18|| , however, that 



upcoming detailed balance experiments, which go through complicated compound nuclear 
states, may provide limits that are better than these by two orders of magnitude. Even 



though recent theoretical work |19| , |20| indicates that one cannot expect g p to be much larger 
than 10~ 8 , it remains worthwhile to translate limits from other experiments into limits on 
g p . Theoretical expectations are easily and often confounded, and it's important to know 
which of the many experiments reported in the literature (and still to come) have the best 
chance of actually seeing time reversal violation. 
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FIGURES 



FIG. 1. The piece of EE2 arising from E2 transitions with Jf = | <— nJj = | using the FPBPN 
force in the large space. The solid line is the sum from eq. (Q). The points correspond to the 
individual points in the sum. The dashed line is the individual B(E2) in e 2 fm 4 divided by a factor 
of 10 5 . It is apparent that ee2 is well converged at low excitation energies. 

FIG. 2. The piece of emi arising from Ml transitions with Jf = | <— nJi = | using the 
FPBPN force in the large space. The solid line is the sum from eq. (Q). The points correspond to 
the individual points in the sum. The dashed line is the individual B(M1) in nuclear magnetons 
divided by a factor of 100. Emi is well converged at excitation energies above 20 MeV. 
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TABLES 



TABLE I. The absolute values of the E2 and Ml matrix elements for the FPVH and FPBPN 
forces compared to the experimentally determined values. The total calculated E2 and Ml strengths 
are also included. 



Transition 


Experiment 


FPVH(small space) 


FPVH(large space) 


FPBPN(large space) 


Ml jhi -s 


■ l)i 

2> 1 


0.113 


0.113 


0.113 


.0.113 


E2 §)i - 


in 


13.35 


7.091 


6.866 


12.88 


Ml §) 2 -> 


' 2)l 


0.344 


1.202 


0.806 


0.858 


E2 §) 2 -> 


3 \ 
2)1 


27.83 


27.62 


36.38 


37.73 


Ml §) 2 -> 


■ 2)1 


0.298 


0.141 


0.018 


0.126 


E2 f ) 2 - 


5 \ 
2)1 


3.208 


11.66 


14.12 


13.87 


B(M1) |)i 


->7l| 




1.154 


11.56 


10.57 


B(E2) §)! • 






211.7 


440.7 


480.5 


B(M1) f)! 


->n| 




0.598 


5.18 


3.63 


B(E2) §)! - 


->n| 




43.8 


97.0 


101.5. 



TABLE II. The phases, with g p = 1 in Eq. (||) and multiplied by a factor of 10 3 , computed with 
the FPVH and FPBPN interactions. The number of Lanczos iterations is listed to illustrate the 
convergence of the phases. 200 iterations were performed only for the cases listed. Each of columns 
3-6 corresponds to one of the terms in Eq. (||); for example, the heading <S£ 2 (n3/2) corresponds to 
the first term of the equation with Jj = | and nJj = |. The last column contains the final phase 
calculated according to Eq. (||). 



Force(Space) 


Lanczos 


££ 2 (n3/2) 


££2(n5/2) 


£Afi(n3/2) 


£jtfl(«5/2) 


£E2 - SMI 




Iterations 












FPVH(small) 


Complete 


-5.08 


-5.15 


7.04 


10.10 


-27.6 


FPBPN (small) 


Complete 


-6.50 


-5.84 


7.94 


13.03 


-33.3 


FPVH(large) 


100 


-3.239 


-2.322 


-0.814 


12.681 


-17.4 


FPVH(large) 


200 


-3.239 




-0.691 






FPBPN (large) 


60 


-2.436 


-2.306 


1.913 


8.5052 


-15.2 


FPBPN (large) 


100 


-2.425 


-2.298 


2.135 


8.4813 


-15.3 


FPBPN(large) 


200 






1.769 







s 



0.0040 



0.0020 - 



CM 
LU 
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0.0000 - 



-0.0020 - 
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E2, nJ, = 5/2 
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E(5/2) n - £(3/2), (MeV) 



20.0 




-0.005 1 1 1 1 

0.0 10.0 20.0 
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